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Проводится в векторно-матричной форме синтез управляющего сигнала нестационарных динамических систем по желаемой ха- 
рактеристике выходного сигнала и предложен алгоритм оценки их параметров методом изображающих векторов. Это оператор- 
ный метод, который всякой временной функции на конечном промежутке времени ставите соответствие п- мерный вектор, а ли- 
нейному оператору ~ матрицу (пхп). Дальнейшие преобразования, необходимые для оценки параметров и управления систе- 
мы, ведутся численными методами. Все это позволяет успешно использовать вычислительную технику, а окончательный резуль - 
тат на основании формулы обращения записывать в аналоговой форме. 


Современные технологические процессы и про- 
изводства характеризуются многофакторностью и 
сложными зависимостями между параметрами. И 
поэтому широкое распространение систем с пере- 
менными параметрами в области автоматического 
управления, в информационно-измерительных си- 
стемах, а также необходимость более глубокого ко- 
личественного и качественного изучения процес- 
сов, протекающих в таких системах, приводит к ин- 
тенсивной разработке цифровых методов синтеза и 
анализа подобных систем и объектов. 

Качество управления объектами в динамике 
определяется многими факторами. Это наличие 
возмущений и их характер, тип объекта - стацио- 
нарный, нестационарный, линейный, нелинейный. 
Тем не менее, доминирующим фактором, во многом 
определяющим эффективность решения задачи 
управления, является положенная в основу исследо- 
вания математическая модель исследуемого объекта 
или процесса. Для решения задачи управления и 
идентификации нестационарными объектами ис- 
пользован метод изображающих векторов [1—3], ко- 
торый сочетает в себе как цифровые, так и аналити- 
ческие приемы решения поставленной задачи. 

Суть метода изображающих векторов состоит в 
том, что каждой функции//) ставится в однознач- 
ное соответствие вектор Р={{\/ Ъ ... х [ р ). Для функции 
//), определенной на промежутке времени [0, /„], 
имеет место разложение 

до = хлад, (о 

уЫ О 

где / - коэффициенты Фурье; ГД г) — ортонорми- 
рованные смещенные полиномы Чебышева І-го ро- 
да; г=/// 0 - безразмерная независимая переменная. 


Приведем некоторые свойства метода изобра- 
жающих векторов. Операции интегрирования 
функции/ г) соответствует в области изображаю- 
щих векторов умножению ее изображающего век- 
тора на матрицу интегрирования: 

V = ІР+ у (0)е ] /Г 0 (т), (2) 

где I - матрица интегрирования, у(0) - начальные 
условия, е, - единичный вектор, Г 0 (т) - полином 
Чебышева. 

Для многократного интегрирования при нуле- 
вых начальных условиях матрица интегрирования 
возводится в соответствующую степень: 

У = І к Р. 

Аналогичное положение имеет место и для ма- 
трицы дифференцирования В [1]. 

Произведению двух функций 
/г(г) = Дг)/(г) 

в области изображающих векторов соответствует 
соотношение вида 

Н = 2Ц)Р, (3) 

где / - матрица Якоби [1]. 

Таким образом, изображение произведения 
двух функций равно произведению изображений 
матрицы известной функции Дт) на изображаю- 
щий вектор другой. Изображающей матрицей 
условно названа матричная функция 2(1), которая 
получается из заданной функции Дт) заменой ска- 
лярного аргумента г на матрицу /. Ввиду равноз- 
начности двух функций их произведение коммута- 
тивно, т. е. равно произведению матрицы второй 
функции на изображающий вектор первой. Тогда 
выражение (3) запишется как 
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Я = 0 г ёіа§[Дг0, Дт 2 ),...,2(т г ,)]0Д (4) 

где 0 - интерполяционная матрица [1], составлен- 
ная из значений базовых полиномов в узлах интер- 
полирования, р - размерность интерполяционной 
матрицы 0, і(т к ) - значения функции і(т) в нулях 
р-то полинома Чебышева. Для учета интервала раз- 
ложения матрица Якоби /умножается скалярно на 
величину / 0 . Восстанавливается исходная функция 
времени /(г) по изображающему вектору в соответ- 
ствии с формулой обращения 

/(т) = (Г,Т(т)), (5) 

где правая часть имеет смысл скалярного произве- 
дения изображающего вектора на переменный век- 
тор полиномов Чебышева Т(т). 

При описании объекта управления в виде пере- 
даточной функции уравнение, связывающее вы- 
ходной сигнал у(Г) и сигнал управления и(1) в опе- 
раторной форме, запишется как 

У = ЖШ, (6) 

где Щя) - передаточная функция объекта управле- 
ния. 

Изображающий вектор выходного сигнала у (Г) 
на интервале [0, / 0 ] (6) можно записать согласно со- 
отношениям (1,2) как 

ѵ=(( 0 і)Щ( а і)и. (7) 

Разрешая соотношение (7) относительно векто- 
ра управляющего сигнала, получим 

і/=ті 0 іу'(і 0 і)г. (8) 

Для матрицы дифференцирования выражение 
(6) запишется как 

Д = ИД/))У. (9) 

Выходной сигнал у(і) можно задать в виде про- 
извольной функции перемещения объекта упра- 
вления по заданным координатам. Для получения 
изображающего вектора произвольного выходного 
сигнала Тна интервале времени [0, / 0 ] воспользуем- 
ся соотношением перехода от точечного изобра- 
жающего вектора к изображающему вектору У сле- 
дующего вида 

У = 2/рО Т У(т к ), (10) 

где У(т к )={у(т к )} - точечный изображающий век- 
тор, представляющий собой вектор, элементы ко- 
торого есть значения функции у(1) в нулях первого 
из отброшенных полиномов Чебышева. С учетом 
точечного изображающего вектора (10) соотноше- 
ния (8, 9) запишутся как 

II = 7/ р\Ѵ (Г 0 І)'О т У(т к ), 

II = 2/ р(Ѵ (О)0 т У(т к ). 

Рассмотрим динамическую систему, которая 
описывается дифференциальным уравнением вида 

А 0 (і)У (п \() + МОУ (п ~ 1) (0 + - + 

+А 1 „_ І) (1)У(0 + А п (1)У (() = [/, (11) 


где /=0 ,й - переменные коэффициенты, не- 
прерывные на интервале [0, / 0 ]. 

На основании правил перехода к изображаю- 
щим векторам (1-4) дифференциальное уравнение 
(6) запишется следующим образом 

\А 0 (1)0" + А, (І)0"~ 1 +... + (1)0 + А п (/)] = V 

или С[ДУ),/>(У)]Г=С/. (12) 

При расчете реакции динамических нелинейных 
нестационарных систем управления на управляю- 
щее воздействие заданные нелинейные коэффици- 
енты по методу изображающих векторов переводят в 
точечные изображающие векторы. Таким образом, 
имея значения неизвестной функции в определен- 
ных точках, узлах интерполирования, можно прове- 
сти определенную нелинейную операцию над этими 
точками, например: операцию вычисления модуля 
функции; нелинейностей, заданных в виде графика 
насыщения; различных тригонометрических и сте- 
пенных зависимостей в этих точках. Далее получен- 
ный преобразованный точечный изображающий 
вектор переводят в изображающий вектор соответ- 
ствующей нелинейной операции. 

Рассмотрим систему, которая описывается не- 
линейным дифференциальным уравнением сле- 
дующего вида: 

Г(і,у) + Щі, У )} = 1/, (13) 

где левая часть уравнения (13) состоит из линейной 
части Д і,у) и нелинейной Ь(і,у). 

По методу изображающих векторов рассмотрим 
вычисление неизвестной функции у (7) при дей- 
ствии нелинейного оператора /. Переходим к то- 
чечным изображающим векторам, тогда нелиней- 
ная часть соотношения ( 1 3) запишется как 

= 0 г ёіа§[1у(т,)]0, 

где К - оператор перехода от точечных изображаю- 
щих векторов функций, стоящих под знаком опе- 
ратора, к изображающим векторам этих функций. 

В векторно-матричной форме соотношение 
(13) запишется в следующем виде 

0(^,У) + ^ т йт ё [^У к \^ = ^. 

Восстанавливается управляющая функция вре- 
мени 11(1) по изображающему вектору (8, 9, 12, 13) 
и в соответствии с формулой обращения 

Ѵ(т) = (ІІ,Т(т)), (14) 

Таким образом, при описании объектов упра- 
вления различными методами (передаточная 
функция или соответствующее дифференциальное 
уравнение) можно получить вектор управляющего 
сигнала в зависимости его от выходного сигнала У, 
преобразовать по соотношению (14) и подать на 
вход системы сигнал в аналоговой форме. 

Задача оценки параметров динамической си- 
стемы сводится к нахождению: 

• формальной структуры модели; 
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• конструктивных параметров объекта управления; 


• неизвестных параметров модели системы упра- 
вления. 


Динамику системы управления успешно опи- 
сывают с помощью дифференциальных и инте- 
гральных уравнений [4] 


Ув,— — у(і) = 'У'Ь. 

П ' аг" й ' м~ т 


сі і ~ 


и(Г), 


(15) 


Формальная структура системы управления на 
основании многих источников определена в виде 
дифференциального уравнения (15). Оператор 
идентификации объекта управления запишем в ви- 
де передаточной функции 

Ж(*) = 4г1 (16) 

А(а) 


Таким образом, задача оценки параметров пере- 
даточной функции по виду переходного процесса 
системы управления Н{і) сводится к оценке кон- 
структивных параметров передаточной функции 
(16), при этом т<п, и нахож дению неиз вестных па- 
раметров модели а,(/ = 0,« — 1) и Ь]Ц=0,т). 

Для решения поставленной задачи воспользу- 
емся методом изображающих векторов, для этого 
переходной процесс к(і) преобразуем в изображаю- 
щий вектор по соотношениям (1,2): 

}• (17) 

По вектору (17) и из формулы обращения (5) 
получаем аналитическую зависимость переходного 
процесса в виде полинома порядка р: 

/7(0 = (Я, тт (18) 

По аналитической зависимости (18) функции 
к(() находим время переходного процесса / 0 от мо- 
мента включения системы до момента, когда мо- 
дуль отклонения переходного процесса от устано- 
вившегося значения не превосходит заданной ве- 
личины зоны нечувствительности. 

Для нахождения весовой функции системы \ѵ(1) 
воспользуемся дифференциальной связью между 
весовой и переходной функциями 

(і9) 


В области изображающих векторов соотноше- 
ние (19) запишется как 

ѵг = ОН. 


По соотношению (5) получили аналитический 
вид весовой функции 

*К0 = (> ѵ,Т({)). 


Для оценки нулей передаточной функции вво- 
дится постоянный коэффициент с> 1 и на интерва- 
ле от 0 до / 0 по соотношению и'(/)іѵ(с/)<0 оценивает- 
ся порядок числителя передаточной функции (16), 
т. е. количество подынтервалов, на которых нахо- 
дятся нули весовой функции и»(/). 


По весовой функции (19) получаем числовую ха- 
рактеристику передаточной функции, для этого вос- 
пользуемся интегральным преобразованием Лапласа 

со 

^(&) = |н</)ехр(-аОЛ, (20) 

0 


где §і - вещественный параметр на интервале [0,4]; 
Щ&)- оператор системы управления. 

Оператор системы управления (15) для веще- 
ственной переменной § запишется как 




ь т 8'"+ь т _х- 1 +.~+ь 1§+ ь 0 

8 п +а„_ 1 8 п - 1+ - +а і8 + а о 


( 21 ) 


Для оценки конструктивных параметров п, т 
воспользуемся уравнениями (20, 21) и предельным 
соотношением 


ц т ЖЙ1 = с "- т 

*-к. РѴ(С8) 


( 22 ) 


Из полученного соотношения (22) находится 
оценка конструктивных параметров 

Іпс"-'" 

п - т = . 

Іпс 

В результате введенных допущений при оценке 
параметров п, т получается вещественное число, со- 
держащее целую часть и мантиссу. Мантиссу полу- 
ченного выражения принимаем за единицу и при- 
бавляем к целой части. Таким образом, задача на- 
хождения конструктивных параметров и, т решена. 

Для нахождения коэффициентов передаточной 
функции необходимо решить систему линейных 
алгебраических уравнений порядка л+м+1 следую- 
щего вида: 

д -1 т 

ДтьНЩьѢаЖ = о, 

/-0 у - 0 

к = 0,и + я7 + 1. 

По найденным коэффициентам а ] ,Ь і и кон- 
структивным параметрам пит записываем переда- 
точную функцию исследуемой системы. 

Необходимо отметить, что переходной процесс 
(17) может быть на выходе нестационарной и нели- 
нейной динамических систем. Изложенный метод 
оценки параметров динамических систем позволя- 
ет аппроксимировать эти системы передаточными 
функциями линейных динамических систем. 

Предложенный метод моделирования неста- 
ционарных систем управления позволяет получить 
сигнал управления в зависимости от выходного 
сигнала У и для исследуемой (идентифицируемой) 
системы управления. Действительно, по получен- 
ной передаточной функции идентифицируемой 
системы по соотношениям (8) или (9) и (14) полу- 
чаем управляющий сигнал Щт), который и можно 
подать на вход исследуемой системы управления. 


46 




Естественные науки 


СПИСОК ЛИТЕРАТУРЫ 

1. Осипов В. М. Основы метода изображающих векторов и линей- 
ное преобразование сигналов // В сб.: Вопросы программиро- 
вания и автоматизации проектирования. - Вып. 1. - Томск: 
Изд-во Томского ун-та, 1971. - С. 1-13. 

2. Осипов В.М., Шалаев Ю.Н. Решение линейных дифферен- 
циальных уравнений с переменными коэффициентами на 


АВМ методом изображающих векторов // Известия вузов. 
Приборостроение. - 1977. - № 12. - С. 43-47. 

3. Шалаев Ю.Н. Применение метода изображающих векторов к 
решению краевых задач для обыкновенных дифференциаль- 
ных уравнений // В сб.: Автоматизация управления и АСУ ТП. 
- Томск: Изд-во ТПУ, 1977. - С. 101-105. 

4. Бесекерский В.А., Попов Е.П. Теория систем автоматического 
регулирования. - М.: Наука, 1972. - 768 с. 


ѵпк чта п пі 




